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We study the dynamics of the spiral magnet under the charge current by solving the Landau- 
Lifshitz- Gilbert equation numerically. In the steady state, the current j induces (i) the parallel shift 
of the spiral pattern with velocity v — (f3/a)j (a, f3: the Gilbert damping coefficients), (ii) the 
uniform magnetization M parallel or anti-parallel to the current depending on the chirality of the 
spiral and the ratio (3 /a, and (iii) the change in the wavenumber k of the spiral. These are ana- 
lyzed by the continuum effective theory using the scaling argument, and the various nonequilibrium 
phenomena such as the chaotic behavior and current-induced annealing are also discussed. 
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The current-induced dynamics of the magnetic struc- 
ture is attracting intensive interests from the viewpoint of 
the spintronics. A representative example is the current- 
driven motion of the magnetic domain wall (DW) in fer- 
romagnets [l|, |2|. This phenomenon can be understood 
from the conservation of the spin angular momentum, 
i.e., spin torque transfer mechanism [sl, [j, IS, [y, [fl]. The 
memory devices using this current-induced magnetic DW 
motion is now seriously considered [8\ . Another example 
is the motion of the vortex structure on the disk of a 
ferromagnet, where the circulating motion of the vortex 
core is sometimes accompanied with the inversion of the 
magnetization at the core perpendicular to the disk [9]. 
Therefore, the dynamics of the magnetic structure in- 
duced by the current is an important and fundamental 
issue universal in the metallic magnetic systems. On the 
other hand, there are several metallic spiral magnets with 
the frustrated exchange interactions such as Ho metal 
[lO|, [ll], and with the Dzyaloshinskii-Moriya(DM) inter- 
action such as MnSi [H, [ll, [13 , (Fe,Co)Si [15], and FeGe 
[.16] . The quantum disordering under pressure or the non- 
trivial magnetic textures have been discussed for the lat- 
ter class of materials. An important feature is that the 
direction of the wavevector is one of the degrees of free- 
dom in addition to the phase of the screw spins. Also the 
non-collinear nature of the spin configuration make it an 
interesting arena for the study of Berry phase effect [17], 
which appears most clearly in the coupling to the current. 
However, the studies on the current-induced dynamics of 
the magnetic structures with finite wavenumber, e.g., an- 
tiferromagnet and spiral magnet, are rather limited com- 
pared with those on the ferromagnetic materials. One 
reason is that the observation of the magnetic DW has 
been difficult in the case of antiferromagnets or spiral 
magnets. Recently, the direct space- time observation of 
the spiral structure by Lorentz microscope becomes pos- 
sible for the DM induced spiral magnets [15 , 16] since the 
wavelength of the spiral is long (~ lOOnm). Therefore, 
the current-induced dynamics of spiral magnets is now 
an interesting problem of experimental relevance. 



In this paper, we study the current-induced dynamics 
of the spiral magnet with the DM interaction as an ex- 
plicit example. One may consider that the spiral magnet 
can be regarded as the periodic array of the DW's in fer- 
romagnet, but it has many nontrivial features unexpected 
from this naive picture as shown below. 

The Hamiltonian we consider is given by [l8| 
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where J > is the exchange coupling constant and 7 is 
the strength of the DM interaction. The ground state of 
H is realized when S{f) is a proper screw state such that 



S{r) = S{ni cos k - f-\- n2 sin /c • r). 



(2) 



where the wavenumber k = usIjI/J, and ni{i = 1,2,3) 
form the orthonormal vector sets. The ground state en- 
ergy is given by —VS'^'y'^ /2J where V is the volume of 
the system. The sign of 7 is equal to that of (ni x 712) -^3, 
determining the chirality of the spiral. 

The equation of motion of the spin under the current 
is written as 

'S = ^B,sxS-^{j-V)S+^0x{j-V)S+^Sx's 

where ^eff = —SH/dS is the effective magnetic field and 
a, P are the Gilbert damping constants introduced phe- 
nomenologically [19, 20]. 

We discretize the Hamiltonian Eq.(l) and the equation 
of motion Eq.(3) by putting spins on the chain or the 
square lattice with the lattice constant a, and replacing 
the derivative by the difference. The length of the spin 
\Si\ is a constant of motion at each site z, and we can 

easily derive H = ^ - S = -a\S\^ from Eq.(3), i.e., the 
energy continues to decrease as the time evolution. 

We start with the one-dimensional case along x-axis 
as shown in Fig.l. The discretization means replacing 

dxS{x) by (6i+i -6i_i)/2a, and dlS{x) by (6^+1 -26^ + 



Si-i)/o?. We note that the wavenumber which mini- 
mizes Eq.([T]) is /c = /co = arcsin(7/J) on the discretized 
one-dimensional lattice. Numerical study of Eq.(3) have 
been done with q/j^b/^ = 1, 2e = 1, 6* = 1, a = 1 
J = 2, and 7 = 1.2. In this condition, the wavelength 
of the spiral A = 27r//co ~ 11.6 is long compared with 
the lattice constant a = 1, and we choose the time scale 
At/(1 + a^) = 10~^. We have confirmed that the re- 
sults do not depend on At even if it is reduced by the 
factor 10~^ or 10~^. The sample size L is 10"^ with the 
open boundary condition. As we will show later, the 
typical value of the current is j ~ 27 and in the real 
situation with the wavelength A[nm], the exchange cou- 
pling constant J[eV] and the lattice constant a[nm], it is 
j ^ 3.2 X 10^V/(Aa)[A/m2]. Substituting J = 0.02, 
A = 100, a = 0.5 into above estimatation, the typi- 
cal current is 10^^[A/m^], and the unit of the time is 
At = J/h^30[ps]. 

The Gilbert damping coefficients a, /3 are typically 
10~^ ~ 10~^ in the realistic systems. In most of the cal- 
culations, however, we take a = 5.0 to accelarate the con- 
vergence to the steady state. The obtained steady state 
depends only on the ratio (3/ a except the spin configura- 
tions near the boundaries as confirmed by the simlations 
with a = 0.1. We employ the two types of initial condi- 
tion, i.e., the ideal proper screw state with the wavenum- 
ber /co, and the random spin configurations. The differ- 
ence of the dynamics in these two cases are limited only 
in the early stage {t < 5000 At). 

Now we consider the steady state with the constant 
velocity for the shift of the spiral pattern obtained after 
the time of the order of 10^ At. One important issue here 
is the current-dependence of the velocity, which has been 
discussed intensively for the DW motion in ferromagnets. 
In the latter case, there appears the intrinsic pinning in 
the case of (3 = [J] , while the highly nonlinear behavior 
for P/a / [20[. In the special case of /3 = a, the trivial 
solution corresponding to the parallel shift of the ground 
state configuration of Eq.(l) with the velocity v = j is 
considered to be realized [5]. Figure 2 shows the results 
for the velocity, the induced uniform magnetization Sx 
along X-axis, and the wavevector k of the spiral in the 
steady state. The current-dependence of the velocity for 
the cases of /3 = 0.1,0.5Qf, a and 2a is shown in Fig. 
2(a). Figure 2(b) shows the /^/a-dependence of the ve- 
locity for the fixed current j = 1.2. It is seen that the 
velocity is almost proportional to both the current j and 
the ratio P/a. Therefore, we conclude that the velocity 
V = {P/a)j without nonlinear behavior up to the current 
j ~ 27, which is in sharp contrast to the case of the DW 
motion in ferromagnets. The unit of the velocity is given 
by a/ At, which is of the order of 20[m/s] for a ~ 5 [A] 
and At ^ 30 [ps]. In Fig. 2(c) shown the wavevector k 
of the spiral under the current j = 1.2 for various values 
of p/a. It shows a non-monotonous behavior with the 
maximum at (3/ a ^0.2, and is always smaller than the 




FIG. 1: Spin configurations in the spiral magnet (a) in 
the equilibrium state, and (b) under the current. Under 
the current j, the uniform magnetization Sx along the spi- 
ral axis/current direction is induced together with the ro- 
tation of the spin, i.e., the parallel shift of the spiral pat- 
tern with the velocity v. Note that the magnetization 
is anti-parallel/parallel to the current direction with posi- 
tive/negative 7 for /3 < a, while it is reversed for /3 > a, 
and the wavenumber k changes from the equilibrium value. 
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FIG. 2: For the case of 7 = 1.2, the numerical result is shown. 
(a) The steady state velocity v as a function of the current 
j with the fixed values of /3 = 0.1, 0.5a, a, and 2a. (b) The 
velocity v as a function of P/a for a fixed value of the current 
j = 1.2. V is almost proportional to P/a. (c) The wavenumber 
k as a function of p/a. The dotted line shows ko in the 
equilibrium, (d) The uniform magnetization Sx along the 
current direction as a function of P/a for a fixed value of 
j = 1.2. 



wavenumber ko in the equilibrium shown in the dotted 
line. Namely, the period of the spiral is elongated by the 
current. As shown in Fig. 2(d), there appears the uni- 
form magnetization Sx along the x-direction. Sx is zero 
and changes the sign at p/a = 1. With the positive 7 
(as in the case of Fig. 2(d)), Sx is anti-parallel to the 
current j/ /x with P < a and changes its direction for 
P > a. For the negative 7, the sign of Sx is reversed. As 
for the velocity v, on the other hand, it is always parallel 
to the current j. 



Now we present the analysis of the above results in 
terms of the continuum theory and a scaling argument. 
For one-dimensional case, the modified LLG Eq.(3) can 
be recast in the following form: 

S = -JSxdlS-{2-fS^^j)d^S^Sx{aS^Pjd^S). (4) 

It is convenient to introduce a moving coordinates ^(x, t), 
7), and C(^5^) (see Fig(T]) [21]. They are explicitly defined 
through f , y and z as 

C{x^t) = cos{k{x — vt) -\- (j))y -\-sm{k{x — vt) -\- (j))z, 
^{x^t) = —sm{k{x — vt)-\-(j))y-\-cos{k{x — vt)-\-(j))z^ 

and fj = X. We restrict ourselves to the following ansatz: 

S{x, t) = S^fj + ^r^^C(^, t), (5) 

where Sx is assumed to be constant. 

By substituting Eq.([5j) into Eq.(|4j), we obtain v as 



(a) 



a 



(6) 



by requiring that there is no force along fj- and ("- 
directions acting on each spin. In contrast to the DW 
motion in ferromagnets, the velocity v becomes zero when 
/3 ^ even for large value of the current. The numerical 
results in Fig. 2(a), (b) show good agreement with this 
prediction Eq.(6). 

On the other hand, the magnetization Sx along x-axis 
is given by 



^x — 



2j-Jk 



J^ 



(7) 



once the wavevector k is known. Here we note that the 
above solution is degenerate with respect to /c, which 
needs to be determined by the numerical solution. From 
the dimensional analysis, the spiral wave number k is 
given by the scaling form, k = kog{j /{2j)^ /3/a) with 
the dimensionless function g{x^y) and also is Sx through 
Eq.(7). 

Motivated by the analysis above, we study the 7- 
dependence of the steady state properties. In Fig. 3, 
we show the numerical results for k/ko and Sx as the 
functions of j/27 in the cases of /3/a = 0.1,0.5 and 
2. Roughly speaking, the degeneracies of the data are 
obtained approximately for each color points (the same 
P/a value) with different 7 values. The deviation from 
the scaling behavior is due to the discrete nature of the 
lattice model, which is relevant to the realistic situation. 
For /3/a = 0.1 (black points in Fig. 3), k remains constant 
and Sx is induced almost proportional to the current up 
jf/27 ^ 0.4, where the abrupt change of k occurs. For 
P/a = 0.5 (blue points) and P/a = 2.0 (red points), the 
changes of k and Sx are more smooth. A remarkable 
result is that the spin Si on the lattice point i is well 
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FIG. 3: The scaling plot for (a) k/ko where ko is the wavenum- 
ber in the equilibrium without the current, and (b) Sx in the 
steady state as the function of j/2j. The black, blue, and 
red color points correspond to (3/a = 0.1, 0.5 and 2.0, respec- 
tively. The curves in (b) indicate Eq.(7) calculated from the 
k values in (a), showing the good agreement with the data 
points. 



described by Eq.(5) at x = x^, and hence the relation 
Eq.(7) is well satisfied as shown by the curves in Fig. 
3(b), even though the scaling relation is violated to some 
extent. For larger values of j beyond the data points, 
i.e., j727 > 0.75 for p/a = 0.1, ^727 > 1.5 for p/a = 0.5 
and j/27 > 0.9 for P/a = 2.0, the spin configuration is 
disordered from harmonic spiral characterized by a sin- 
gle wavenumber k. The spins are the chaotic funtion of 
both space and time in this state analogous to the turbu- 
lance. This instability is triggered by the saturated spin 
Sx = ±1, occur ing near the edge of the sample. 

Next, we turn to the simulations on the two- 
dimensional square lattice in the xy-plane. In this case, 
the direction of the spiral wavevector becomes another 
important variable because the degeneracy of the ground 
state energy occurs. 

Starting with the random spin configuration, we sim- 
ulate the time evolution of the system without and with 
the current as shown in Fig. 4. Calculation has been done 
with the same parameters as in the one-dimensional case 
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FIG. 4: The time evolution of the z component Sz of the 
spin from the random initial configuration of the 10^ x 10^ 
section in the middle of the sample is shown in the case of 
(a) j = 0, (b) j = 0.3 along the x-axis, (c) j = 0.3 along 
the (l,l)-direction. From the left, t = 10^ At, 1.7 x 10^ At, 
5 X 10^ At. The rightmost panels show the spectral intensity 
\Sz{k, 5 X lO^At)!^ from the whole sample of the size 2^° x 2^° 
in the momentum space k = (kx^ky). 



where 7 = 1.2, /3 = 0, and the system size is 2^^ x 2^^. 

In the absence of the current, the relaxation of the 
spins into the spiral state is very slow, and many dislo- 
cations remain even after a long time. Correspondingly, 
the energy does not decrease to the ground state value 
but approaches to the higher value with the power-law 
like long-time tail. The momentum-resolved intensity is 
circularly distributed with the broad width as shown in 
Fig. 4(a) corresponding to the disordered direction of 
k. This glassy behavior is distinct from the relaxation 
dynamics of the ferromagnet where the large domain for- 
mation occurs even though the DW's remain. Now we 
put the current along the x (Fig. 4(b)) and {x-\-y) (Fig. 
4(c)) directions. It is seen that the direction of k is con- 
trolled by the current also with the radial distribution 
in the momentum space being narrower than that in the 
absence of j (Fig. 4(a)). This result suggests that the 
current j with the density ~ 10^^[A/m^] of the time du- 
ration ~ 0.1[/isec] can anneal the directional disorder of 
the spiral magnet. After the alignment of k is achieved, 
the simulations on the one-dimensional model described 
above are relevant to the long-time behavior. 



To summarize, we have studied the dynamics of the 
spiral magnet with DM interaction under the current j 
by solving the Landau-Lifshitz- Gilbert equation numeri- 
cally. In the steady state under the charge current j, the 
velocity v is given by {/3/a)j (a, /3: the Gilbert-damping 
coefficients), the uniform magnetization is induced par- 
allel or anti-parallel to the current direction, and period 
of the spiral is elongated. The annealing effect especially 
on the direction of the spiral wavevector k is also demon- 
strated. 
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